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Abstract 

Novel  characterization  of  the  thermal  properties  of  batteries  have  been  introduced  by  defining  their  frequency-dependent  thermal 
impedance  functions.  The  thermal  impedance  function  can  be  approximated  as  a  thermal  impedance  spectrum  by  analyzing  the  experimental 
temperature  transient  which  is  related  to  the  thermal  impedance  function  through  Laplace  transformation. 

In  order  to  obtain  the  temperature  transient,  a  process  has  been  devised  to  generate  an  external  heat  pulse  with  heating  wire  and  to  measure 
the  response  of  the  battery.  This  process  is  used  to  study  several  commercial  Li-ion  batteries  of  cylindrical  type.  Thermal  impedance 
measurements  have  been  performed  using  a  potentiostat/galvanostat  controlled  by  a  digital  signal  processor,  which  is  more  commonly 
available  than  a  flow-meter  usually  applied  for  thermal  property  measurements. 

Thermal  impedance  spectra  obtained  for  batteries  produced  by  different  manufacturers  are  found  to  differ  considerably.  Comparison  of 
spectra  at  different  states-of-charge  indicates  an  independence  of  the  thermal  impedance  on  the  charge  state  of  the  battery.  It  is  shown  that  the 
thermal  impedance  spectrum  can  be  used  to  obtain  simultaneously  the  thermal  capacity  and  the  thermal  conductivity  of  the  battery  by  non¬ 
linear  complex  least-squares  fit  of  the  spectrum  to  the  thermal  impedance  model.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

As  batteries  have  been  principally  employed  in  low-power 
applications,  their  thermal  properties  have  not  been  relevant 
to  the  application  design  and  therefore,  have  not  been  the 
subject  of  detailed  research.  This  situation  has  changed 
markedly  with  the  onset  of  new,  high-power  applications 
of  batteries,  such  as  portable  computers,  video-cameras, 
power-tools  and  hybrid  electric  vehicles.  As  the  charge- 
discharge  rates  used  in  such  devices  can  reach  the  15  C  rate 
(full  discharge  within  4  min),  substantial  heat  development 
can  occur  in  a  short  time.  Thus,  advanced  heat  management 
of  batteries  becomes  important,  particularly  for  large  bat¬ 
teries  or  battery  stacks,  to  prevent  overheating  damage  or 
thermal  runaway  [1].  This  is  especially  the  care  with  Li-ion 
batteries  which  are  capable  of  explosion  and  fire  if  over¬ 
heated.  The  requirements  of  the  engineering  of  high-power 
battery  systems  have  resulted  in  increased  attention  being 
paid  to  the  modeling  of  the  thermal  [2-4]  and  the  combined 
electrochemical-thermal  [5,6]  behavior  of  batteries.  For 
such  modeling,  it  is  necessary  to  know  the  experimental 
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values  of  the  thermal  conductivity  and  the  thermal  capacity 
of  each  particular  battery  type.  The  thermal  properties  of 
many  battery  components  have  been  studied  by  means  of  a 
guarded  heat  flow-meter  used  by  Song  et  al.  [7,8],  By  way  of 
example,  theoretical  calculations  of  the  thermal  conductivity 
and  the  capacity  of  an  electrode  stack  of  a  Li-polymer 
battery  are  given  in  [4].  Such  calculations  are  only  possible 
if  the  composition  of  the  battery  stack  is  known  exactly. 
Often,  this  information  is  available  only  to  battery  manu¬ 
facturers.  A  further  problem  arises  from  the  fact  that  the 
thermal  characteristics  of  the  stack  components  can  change 
after  exposure  to  electrolyte  and  the  formation  cycle.  In  situ 
measurements  of  these  characteristics  would  give  more 
precise  information.  Several  recent  papers  have  provided 
experimental  data  on  the  thermal  capacity  of  the  cells, 
measured  by  calorimetric  technique  [9,10].  By  contrast, 
no  in  situ  measured  experimental  data  for  the  thermal 
conductivity  of  Li-ion  batteries  has  been  reported  to  date. 

In  this  study,  we  present  a  new  method  to  measure 
simultaneously  both  the  thermal  conductivity  and  the  capa¬ 
city  of  a  battery  which  employs  commonly  available  elec¬ 
trochemical  equipment  without  the  use  of  a  calorimeter.  The 
new  method  adopts  the  concept  of  the  thermal  impedance  as 
a  complete  characteristic  of  a  system’s  thermal  properties. 
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Such  an  approach  has  been  employed  before  in  different 
areas  [11-13],  The  method  of  measuring  the  thermal  impe¬ 
dance  using  Laplace  transformation  of  the  temperature 
response  transient  to  an  external  or  an  internal  heat  pulse 
is  described  in  detail.  To  justify  the  new  method  experimen¬ 
tally,  the  thermal  impedance  spectra  of  cylindrical  Li-ion 
batteries  (18650  type)  made  by  different  manufacturers 
have  been  obtained.  Especially,  the  dependence  of  thermal 
impedance  on  the  state-of-charge  is  investigated  for  a  Pana¬ 
sonic  Li-ion  battery.  The  specific  thermal  capacity  and  the 
conductivity  are  obtained  for  all  batteries  using  a  non-linear 
fit  of  the  thermal  impedance  spectra  to  a  model  impedance 
function  in  the  frequency-domain,  developed  for  the  case  of 
cylindrical  rod  with  interfacial  heat  development. 


2.  Derivation  of  thermal  impedance  functions 

The  thermal  capacity  and  the  thermal  conductivity  can  be 
obtained  by  direct  analysis  of  the  temperature  transient, 
which  is  measured  after  applying  a  heat  pulse.  Such  analysis 
can  be  performed,  for  example,  by  fitting  time-domain 
experimental  data  to  a  model  function.  Such  a  fit  is  often 
not  practical  where  system  is  so  complicated  that  its 
response  cannot  be  described  analytically  and  has  to  be 
calculated  numerically  using  considerable  computing 
resources.  When  some  unaccounted  thermal  time-constants 
are  present  in  the  system,  the  fit  to  a  beforehand-defined 
model  function  will  not  give  indication  about  their  presence 
or  their  nature.  An  alternative  approach  to  analyze  circuits 
with  unknown  structure  is  adapted  in  electronics.  Repre¬ 
sentation  of  the  transfer  function  of  a  circuit  as  a  frequency- 
domain  spectrum  allows  the  parameters  for  all  known 
elements  of  the  circuit  to  be  obtained  while  detecting  all 
unexpected  elements  in  case  they  are  present. 

The  complex  transfer  function  F(s,  x)  defined  in  a 
Laplace-domain  is  commonly  derived  in  an  intermediate 
step  for  the  temperature  transient  calculation  by  the  Laplace 
method.  If  a  small  temperature  change  is  expected  transients 
can  be  derived  as  T{t )  =  L_1{T’(^,x)},  where  L~l  is  the 
inverse  Laplace  transform,  s  the  complex  frequency  2nf 
i  +  b,  f  the  frequency,  i  the  imaginary  multiplier,  b  the 
arbitrary  real  number,  and  x  the  space  coordinate.  Such 
functions  have  been  derived  for  large  variety  of  shapes  of 
the  heat  conductors  and  heat  application  boundaries  as 
solutions  of  subsidiary  equations  [14]. 

Exploiting  the  mathematical  analogy  between  the  form¬ 
alisms  which  describe  linear  time-invariant  electrical  net¬ 
works  and  thermal  conductors,  we  can  represent  F(s,  .r)  as  a 
product  of  two  functions  I(s)  and  Zt(s,  x).  Here,  the  first 
function  is  the  Laplace  transform  of  heat  excitation  and  the 
second  is  a  complete  description  of  the  thermal  properties  of 
the  system  which  are  dependent  on  the  thermal  conductivity, 
thermal  capacity  and  shape  of  the  thermal  conductor. 

When  the  function  Zt(s,  x )  is  known,  in  a  special  case 
when  b  =  0  and  space  coordinate  is  fixed,  Zt(2nfi)  is 


mathematically  analogous  to  electrical  impedance  and  is 
called  here  the  ‘thermal  impedance  function’.  In  contrast  to 
time-domain  response  functions,  thermal  impedance  func¬ 
tions  have,  in  most  cases,  analytical  expressions.  These 
functions  can  be  evaluated  for  a  set  of  given  frequencies 
fi  as  a  set  of  complex  numbers  Zf,  which  represents  a  thermal 
impedance  spectrum.  The  thermal  impedance  spectrum  can 
be  visualized  as  an  electrical  impedance  spectrum  in  a 
complex  plot.  The  spectrum  in  such  presentation  gives 
indications  about  the  contributions  of  different  parts  of 
the  system  to  heat  transfer  and  allows  detection  of  the 
relative  changes  of  thermal  properties  even  by  simple  visual 
observation.  Thermal  resistance  and  thermal  capacitance 
can  be  geometrically  deduced  from  the  spectrum  in  a 
complex  plot.  For  example,  the  value  of  a  serial  resistor 
can  be  obtained  from  the  highest  frequency  point  of  the 
spectrum,  and  the  value  of  parallel  resistors  from  the  dia¬ 
meter  of  the  semicircle. 

Due  to  the  mathematical  similarity  to  the  electrical  model, 
thermal  impedance  functions  can  be  derived  similarly  to 
electric  impedance  functions,  essentially  using  the  Ohm 
and  Kirchhoff  laws.  In  this  approach,  heat  flow  is  treated 
as  a  current  equivalent  and  temperature  gradient  as  a  voltage 
equivalent.  Accordingly,  the  thermal  resistance  can  be  defined 
through  the  traditional  specific  thermal  conductivity  k  as 
Rt  =  1/kA  (Ks/J),  and  the  thermal  capacitance  through 
the  volume-specific  heat  Cr  =  cAl  (J/K).  Here,  A  is  the 
cross-section  and  /  is  the  length  of  the  elementary  heat 
conduction  region.  By  analogy  with  corresponding  electric 
components,  the  complex  impedance  of  these  components  are 
Z{Rj)  =  Rj  and  Z(CT.  s)  =  1  jCjS  where  .v  is  the  complex 
frequency  2nfi.  Calculation  of  the  complete  thermal  impe¬ 
dance  function  of  a  thermal  conductor  can  be  performed  by 
setting  up  an  “equivalent  thermal  circuit”  and  by  applying  the 
same  rules  as  for  electric  elements  to  find  the  complete 
impedance  between  two  desired  points  of  the  circuit. 

Let  us  consider  the  case  of  a  thin  plate,  with  thickness  /, 
specific  heat  c  and  exposed  to  air  so  that  the  thermal 
conductivity  per  unit  surface  due  to  radiation  cooling  is  k, 
which  is  being  internally  heated.  The  temperature  response 
is  measured  by  a  thermocouple  on  the  surface.  The  thermal 
equivalent  circuit  will  be  a  thermal  resistor  Rrc  =  1  / kA  in 
parallel  with  a  thermal  capacitor  Cc  =  cl  A.  The  thermal 
resistance  of  the  thermocouple  can  be  represented  as  a  serial 
thermal  resistance  Rser.  Such  a  model  gives  a  good  descrip¬ 
tion  of  thermal  behavior  of  the  battery  casing  and  the  heating 
coil.  The  thermal  impedance  of  this  model  is  given  by: 


Z(s) 


1 

—  7 — / - x - b 

(l/Rrc)  +  CcS 


(1) 


Treatment  of  the  function  in  terms  of  thermal  resistance  and 
capacitance  can  be  undertaken  by  impedance  fitting  soft¬ 
ware,  of  the  type  commonly  used  in  electrochemistry,  such 
as  MEISP  developed  by  the  Korea  Kumho  Petrochemical 
Co.  Ltd.,  then  converting  the  obtained  parameters  into 
specific  thermal  values. 
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Fig.  1 .  Thermal  equivalent  circuit  corresponding  to  cylindrical  battery  heated  from  surface  under  radiation  cooling  condition. 


In  many  cases,  an  impedance  function  for  an  electric 
equivalent  circuit  is  known  and  can  be  applied  for  the 
thermal  conduction  model.  For  example,  the  thermal  equiva¬ 
lent  circuit  for  a  set-up,  where  a  constant  heat  flow  is  applied 
to  one  side  of  a  slab  and  another  side  is  thermally  insulated, 
is  equivalent  to  an  open-circuit  terminated  transmission  line 
in  electronics.  Its  impedance  function  is  already  known  as: 


z (s)  =  \j  ^cothi  sf  RtCts)  (2) 

This  equation  (divided  by  a  factor  of  two)  expresses  the 
thermal  impedance  of  a  prismatic  battery  to  which  a  heat 
pulse  is  applied  from  the  two  opposite  wide  sides,  and  where 
the  heat  flow  through  its  thin  sides  can  be  neglected. 

The  thermal  model  for  a  cylindrical  battery,  with  heat 
excitation  applied  from  the  surface  and  with  negligible  heat 
flow  over  the  sides,  is  analogous  to  a  diffusion  model 
developed  for  infinite  cylindrical  rod.  Its  impedance  func¬ 
tion  has  been  derived  as  a  solution  of  the  subsidiary  equation 
in  [14]  and  as  an  electrochemical  impedance  for  diffusion 
towards  the  center  of  the  cylindrical  rod  in  [15]: 


7M  =  Ig(\j2TTCTs )  R 
1  ’  (V  2RTCTs)I\  (V  2RtCts)  T 


(3) 


and  where  70  and  f  are  Bessel  functions  of  the  first  kind,  with 
0  and  1  order,  respectively.  Note  that  transformation  from  a 
cylindrical  coordinate  system  to  a  Cartesian  system  is 
required  to  calculate  specific  thermal  conductivity  and 
capacity  from  the  R  and  C  obtained  from  the  fit.  The 
application  of  the  transformation  results  in  equations  for 
specific  thermal  values,  i.e. 


271 IR 
C 

nr2l 


(4) 

(5) 


where  r  is  the  radius  of  the  battery  and  /  the  length. 

In  addition  to  a  cylindrical  thermal  conductor  which 
represents  the  bulk  of  a  battery  electrode  stack,  the  thermal 
capacity  of  the  battery  casing  together  with  the  heating  coil 
Cc  and  the  thermal  resistance  between  the  casing  and 


electrode  stack  Rc  should  be  included  in  the  thermal  model. 
When  the  battery  surface  is  exposed  to  air,  heat  flow  through 
radiation  and  convection  has  to  be  considered.  Using  the 
thermal  circuit  approach,  this  is  easily  done  by  adding  a 
‘thermal  resistor’  RIC  in  parallel  to  the  impedance  in  Eq.  (3). 
It  can  be  recalculated  give  the  commonly  used  heat  transfer 
coefficient  h  as  h  =  1  / (Rlc2nrl) . 

Additionally,  the  thermal  resistance  between  the  surface 
and  the  thermocouple  should  be  considered  as  a  serial 
thermal  resistance.  The  complete  thermal  equivalent  circuit 
corresponding  to  a  cylindrical  Li-ion  battery  is  shown  in 
Fig.  1.  The  impedance  function  corresponding  to  the  circuit 
is  given  in  Eq.  (6)  as  follows: 


Z(s) 


1 

(lo(V 2RjCjs )/(-\/ 2RtCjs)I\  (\/ 2RjCts))Rt+Rc 


1 


+  Ccs 


+  If 


(6) 


3.  Measurement  of  thermal  impedance  spectrum 


The  thermal  impedance  function  can  be  obtained  directly 
from  time-domain  response  to  heat  excitation  without  any 
presumption  about  internal  structure  of  the  circuit  as  given  in 
Eq.  (7): 


Zt(s ) 


L{T(t)} 

m)} 


(7) 


where  L{  }  denotes  a  Laplace  transform,  T(t)  is  a  tempera¬ 
ture  response  function,  and  i(t)  is  an  excitation  function  in 
the  time-domain.  For  example,  if  i(t)  is  a  heat  pulse,  it  can  be 
expressed  as  heat-flow  amplitude  70  (J/s)  multiplied  by  a 
Heavisde  step-function  (IHt)  defined  as  0  if  t  <  0  and  1 
otherwise,  i(t)  =  70 <P(t).  Its  Laplace  transform  can  be  easily 
calculated  as  L{fd>(t)}  =  Iq/s. 

On  the  other  hand,  the  analytical  form  of  T(t)  is  generally 
not  known.  Information  about  it  is  obtained  from  an  experi¬ 
ment  as  a  set  of  time-dependent  samples  7]  collected  with 
time  interval  5 1  during  sampling  period  At.  The  information 
about  function  Zt(2nfi )  is  valid  in  a  frequency  range  defined 
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by  the  sampling  period  and  interval.  Quantitatively  the 
validity  boundaries  are  given  by  the  Nuquist  theorem, 
/max  =  1  /2 1  and  /m;n  =  1  /At.  A  straightforward  way  to 
calculate  the  impedance  spectrum  based  on  such  set  of 
samples  includes  interpolation  of  the  samples  with  piece- 
wise-linear  or  other  functions,  followed  by  numerical  inte¬ 
gration  of  the  resulting  function  to  evaluate  the  transform 
integral  in  Eq.  (8)  for  s  =  2nfi  +  b  with  frequency /ranging 
from  /max  to  /min,  while  b  is  chosen  so  that  the  integral 
converges.  Different  methods  to  define  b  and  resulting 
approximated  solutions  of  the  transform  integral  have  been 
proposed  in  [16-18]: 

/»oo 

F(s)=  /  f(t)  e_sfdf  (8) 

Jo 

Numerical  integration  and  selection  of  optimal  b  can  be 
avoided  by  using  an  analytical  Laplace  transform  of  a 
piecewise-linear  function.  All  methods  based  on  interpola¬ 
tion  and  direct  use  of  Eq.  (8)  are,  however,  quite  sensitive  to 
noise  which  is  commonly  present  when  7)  are  acquired  in 
experiment,  and  gives  unacceptable  results  at  a  usually 
attainable  noise  level,  as  shown  in  [19].  Alternatively,  con¬ 
version  can  be  accomplished  by  fitting  the  experimental  data 
to  a  physically  relevant  function  which  has  sufficient  degrees 
of  freedom  to  describe  all  possible  system  responses.  The 
parameter  obtained  can  be  substituted  into  an  analytical 
expression  of  the  Laplace  transform  of  such  a  function.  Such 
procedure  called  carrier  function  Laplace  transform  (CFLT) 
is  found  in  [19]  to  exhibit  superior  noise-rejection  capability. 
CFLT  is  therefore  used  here  to  calculate  L  { T( t ) }  which  is 
then  used  to  obtain  the  thermal  impedance  spectrum.  The 
sum  of  an  adaptively  chosen  number  of  exponents  given  in 
Eq.  (9)  is  used  as  a  carrier  function  in  this  case: 

/(0  =  £*.-e-(,/T,)  (9) 

/= 0 


Laplace  transform  is  as  follows: 


n 


ns)  =  £ 

i'=0 


h 

(*+(  1/T.O) 


(10) 


Once  the  function  parameters  k,  and  t ;  are  obtained  from  a  fit 
of  the  experimental  temperature  transient,  the  thermal  impe¬ 
dance  can  be  calculated  combining  Eq.  (10)  with  Eq.  (7),  i.e. 


zw  = 

1=0 


(s+  (1/t,-)) 


(ID 


where  I  (W)  is  the  magnitude  of  the  heat  flow  applied  to  the 
surface  of  the  heat  conductor  under  test  at  time  0,  and  .v  is  2tt.fi 
where  the  frequency/can  be  selected  between  l/8r  and  1/A t. 


4.  Experimental  set-up 

To  demonstrate  the  applicability  of  the  proposed  method 
in  a  real  experimental  environment  and  to  obtain  thermal 


Table  1 

Characteristic  parameters  of  batteries  subjected  to  thermal  impedance 
measurement 


Manufacturer 

Radius/ 
height  (cm) 

Weight  (g) 

Discharge 
capacity  (mAh) 

Moly  Energy 

0.9/6.5 

42.72 

1836 

Panasonic 

0.9/6.5 

39.47 

1536 

Hitachi 

0.9/6.5 

41.57 

1263 

Samsung 

0.9/6.5 

40.93 

1436 

impedance  spectra  for  Li-ion  batteries,  we  used  our  pro¬ 
prietary  potentiostat/galvanostat  with  a  digital  signal  pro¬ 
cessor  (DSP)  for  signal  generation  and  data  acquisition.  The 
system  is  known  as  the  BPS  1000FL  Powergraphy  Battery 
Parameterization  System. 

Cylindrical  Li-ion  batteries  of  the  18650  type  with  a 
specified  capacity  1500  mAh  and  manufactured  by  Sony, 
Panasonic,  LG  and  Samsung  were  obtained  from  commer¬ 
cial  battery  packs.  Detailed  characteristics  of  the  batteries 
are  given  in  Table  1. 

In  order  to  apply  a  constant  heat  flow  to  the  surface  of  the 
batteries,  a  heating  band  made  of  Ni-Fe  alloy  with  a  total 
resistance  of  around  5  Q  and  a  negligible  change  in  this 
value  over  the  experimental  temperature  range  were  tightly 
wound  around  the  battery,  as  shown  in  Fig.  2.  A  very  thin 
contact  wire  was  soldered  on  to  the  heating  band  to  provide 
high  electric  conductivity  but  only  very  low  stray  thermal 
conductivity.  The  opposite  ends  of  the  contact  wire  were 
connected  to  the  positive  and  the  negative  terminals  of  the 
BPS  1000FL,  respectively.  A  thermocouple  was  fixed  at  the 
mid-point  of  the  battery  so  that  it  touched  the  surface  of 
battery  but  not  the  heating  band.  Readings  from  the  thermo¬ 
couple  were  performed  by  the  BP1000FL  and  transmitted 
to  the  control  PC.  The  battery  was  compressed  between  two 
plates  of  thermal  insulator  (porous  polystyrene)  in  order 
to  minimize  heat  conduction  along  the  battery  axis  and 
through  its  sides.  To  increase  the  transfer  of  heat  from  the 
battery  surface  and,  thereby,  decrease  the  largest  thermal 


Fig.  2.  Experimental  set-up  for  thermal  impedance  measurement  on 
cylindrical  battery  using  external  heating  coil  for  generating  heat  pulse. 
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time-constant  of  the  system,  a  2  W  CPU-cooler  fan  was 
placed  1  cm  from  the  battery  to  provide  convection  cooling. 
Experiments  were  performed  in  an  incubator  at  20  °C  with 
an  average  temperature  fluctuation  of  0.05  °C. 

Before  each  measurement,  the  resistance  of  the  coil  was 
measured  at  a  frequency  of  1  kHz  using  the  impedance 
measuring  capabilities  of  the  BPS  1000FL.  With  the  exact 
value  of  resistance,  the  amount  of  current  to  apply  a  heat 
pulse  of  1  W  to  the  surface  of  the  cell  was  calculated.  The 
current  was  applied  to  the  heating  band  with  a  rise  time  of 
40  ps.  The  temperature  response  was  measured  with  a 
sampling  rate  5 1  =  20  ms  during  period  of  At  =  1000  s. 
Calculations  of  the  impedance  spectra  were  performed 
automatically  by  the  BPS  1000FL  control  software  using 
the  algorithm  described  in  Section  1.  Parameters  k,  were 
found  by  linear  fit  to  eight  exponents  with  time-constants  f, 
logarithmically  selected  between  /max  =  l/8f  and 
/min  =  1  j At.  The  fit  was  performed  using  the  normal  equa¬ 
tion  of  the  Gauss-Jordan  elimination  method  [20]. 


5.  Results  and  discussion 

The  thermal  impedance  of  the  heating  coil  used  for  the 
subsequent  experiments  was  measured  between  0. 1  Hz  and 
1  mHz  in  order  to  obtain  the  specific  heat  capacity  of  the  coil 
which  will  be  used  for  correction  of  battery  data.  A  heat 
pulse  of  0.5  W  was  applied  to  the  coil  for  1000  s.  The 
temperature  transient  data  were  converted  to  a  thermal 
impedance  spectrum  as  described.  The  resulting  thermal 
impedance  spectra  together  with  the  fit  obtained  using  Eq.  ( 1 ) 
are  shown  in  Fig.  3.  The  MEISP  non-linear  least-squares 
fitting  software  developed  by  Kumho  Petrochemical  Co. 
Ltd.  was  used  for  analysis  of  the  spectra.  A  good  fit  of 


the  spectrum  by  the  thermal  impedance  function  indicates 
the  correctness  of  the  chosen  thermal  equivalent  circuit 
model.  The  thermal  capacity  of  the  coil  obtained  by  the 
fit  is  4.12  J/K. 

In  order  to  separate  the  influence  of  the  battery  casing 
from  the  thermal  behavior  of  the  electrode  stack,  the  thermal 
impedance  spectrum  of  the  empty  casing  was  measured.  The 
thermal  capacity  Cc  of  the  casing  together  with  the  heating 
coil  is  found,  by  the  fit  of  the  spectrum,  to  be  13.74  J/K. 

The  thermal  impedance  spectrum  of  the  Panasonic  battery 
was  measured  over  the  wide  frequency  range  of  0.1  Hz- 
100  pHz  to  check  the  validity  of  the  proposed  thermal  model 
and  to  obtain  the  thermal  capacity  and  the  thermal  con¬ 
ductivity  of  the  battery.  The  resulting  temperature  transient 
is  shown  in  Fig.  4  and  the  impedance  spectrum  together  with 
the  model  fit  is  given  in  Fig.  5  as  a  complex  plot.  Two 
semicircles  are  visible  on  the  plot.  The  depressed  semicircle 
in  the  high-frequency  part  corresponds  to  the  thermal  capa¬ 
city  of  the  battery  casing  in  parallel  with  thermal  resistance 
between  the  casing  and  the  electrode  stack.  The  semicircle  is 
depressed  by  the  thermal  resistance  of  the  electrode  stack. 
The  low-frequency  semicircle  corresponds  to  the  thermal 
capacity  of  the  electrode  stack  in  parallel  with  the  thermal 
resistance  of  the  interfacial  cooling. 

The  parameters  of  the  thermal  model  were  obtained  by 
fitting  the  spectrum  to  the  thermal  impedance  function  in 
Eq.  (5).  The  casing  capacity  Cc  obtained  has  been  used  as  a 
fixed  parameter  in  the  fit.  The  good  fit  result  indicates  that 
the  proposed  thermal  model  satisfactorily  describes  the 
thermal  behavior  of  the  cylindrical  Li-ion  battery.  All  para¬ 
meters  resulting  from  the  fit  are  listed  in  Table  2.  As 
expected,  thermal  resistance  of  the  thermocouple  introduces 
a  delay  in  the  temperature  onset  and  is  therefore  represented 
as  a  negative  resistance  in  the  thermal  equivalent  circuit.  The 
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Fig.  3.  Thermal  impedance  spectrum  of  heating  coil,  measured  in  frequency  range  0.2  Hz-1  mHz.  Fit  line  corresponds  to  the  thermal  impedance  function  for 
model  described  by  Eq.  (1). 
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Fig.  4.  Temperature  transient  in  response  to  1  W  heat  pulse  applied  to  surface  of  Panasonic  Li-ion  battery. 


Fig.  5.  Thermal  impedance  spectrum  of  Panasonic  Li-ion  battery  in  fully- 
charged  state,  measured  between  0.1  Hz  and  100  pHz.  Fit  line  corresponds 
to  the  thermal  model  described  by  Eq.  (5)  with  parameter  values  given  in 
Table  2. 

thermal  capacity  of  the  metal  casing  is  found  to  be  respon¬ 
sible  for  one-quarter  of  the  total  thermal  capacity.  A  better 
thermal  contact  between  the  casing  and  the  electrode  stack 
can  be  important  for  improving  the  heat  dissipation  from 
battery,  since  the  thermal  resistance  Rc  between  the  casing 


and  the  stack  is  found  to  be  one-third  of  the  total  thermal 
resistance  of  the  battery.  The  thermal  resistance  at  the 
battery  surface  due  to  radiation  and  convection  cooling 
Rrc  is  the  largest  bottleneck  in  battery  heat  dissipation  in 
spite  of  air  circulation  with  a  fan.  Other  contributions  to  total 
thermal  impedance  area,  however,  of  the  same  order  of 
magnitude.  In  the  case  of  active  cooling  with  water,  as 
considered  for  high-power  applications  such  as  electric 
vehicles,  the  thermal  resistance  of  the  electrode  stack  can 
become  dominant. 

Specific  values  of  the  thermal  parameters  have  been 
calculated  using  information  about  battery  dimensions 
and  weight  in  Table  1  and  Eq.  (4).  To  calculate  the  specific 
capacity  of  the  casing,  the  thermal  capacity  of  the  coils 
subtracted,  and  a  casing  weight  of  8.9  g  was  used.  All 
specific  values  are  given  in  Table  3.  The  specific  thermal 
capacity  of  the  electrode  stack  is  1.9  J/(K  g)  and  is  notice¬ 
ably  higher  than  that  the  battery  previously  reported  as 
1.07  J/(K  g)  in  [9].  The  thermal  capacity  of  the  casing  is 
not  considered  separately  in  this  paper.  The  thermal  con¬ 
ductivity  of  the  electrode  stack  perpendicular  to  the  elec¬ 
trode  surface  is  found  to  be  1 .4  W/(K  m),  which  is  almost 
five  times  larger  than  that  for  a  Li-polymer  battery  calculated 
in  [4]  on  the  basis  of  the  known  thermal  conductivity  of  the 
component  materials.  This  considerable  difference  probably 


Table  2 

Thermal  equivalent  circuit  parameters  obtained  by  fitting  thermal  impedance  spectra  of  Li-ion  batteries  (Fig.  6)  to  thermal  model  in  Fig.  1  described  by 
Eq.  (5) 


Manufacturer 

RSc,  (K/W) 

Rrc  (K/W) 

Cc  (J/K) 

Rc  (K/W) 

CT  (J/K) 

Rt  (K/W) 

Moly  Energy 

-0.05 

2.49 

13.47 

0 

61.5 

2.80 

Panasonic 

-0.1 

1.98 

13.74 

0.90 

58.4 

1.78 

Hitachi 

-0.05 

3.23 

13.74 

1.01 

37.68 

0.86 

Samsung 

-0.02 

1.94 

13.74 

0 

76.73 

2.07 
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Table  3 

Heat  transfer  coefficients,  specific  heat  and  thermal  conductivity  of  cylindrical  Li-ion  batteries,  obtained  by  fit  of  their  thermal  impedance  spectra  to  thermal 
impedance  function  in  Eq.  (5)  corresponding  to  thermal  model  given  Fig.  1 


Manufacturer 

h,c  (W/(K  m2)) 

Cc  (J/(K  g)) 

hc  (W/(K  m2)) 

Ct  (J/(K  g)) 

kt  (W/(K  m)) 

Moly  Energy 

109.3 

1.1 

- 

1.8 

0.9 

Panasonic 

137.1 

1.1 

302.1 

1.9 

1.4 

Hitachi 

84.5 

1.1 

269.3 

1.2 

2.8 

Samsung 

140.2 

1.1 

- 

2.4 

1.2 

Specific  values  are  obtained  using  conversion  equations  given  in  Eq.  (4). 


arises,  apart  from  the  use  of  liquid  electrolyte,  from  the  fact 
that  a  wound  electrode  stack  is  used  in  commercial  batteries 
(both  cylindrical  and  prismatic)  and  not  the  parallel  elec¬ 
trode  arrangement  adopted  in  the  calculations.  A  wound 
electrode  stack  design  results  in  a  thermal  conduction  path¬ 
ways  in  radial  direction  through  highly  thermally  conductive 
current-collectors.  The  presence  of  such  unaccounted  ther¬ 
mal  pathways  can  explain  a  much  higher  average  thermal 
conductivity,  as  found  in  experiment.  The  value  obtained  for 
the  transfer  coefficient,  viz.  137W/(Km2)  is  rather  high 
compared  with  25  W/(K  m2)  which  is  usually  assumed  for 
radiation  and  active  convection  cooling  [4],  This  is  probably 
due  to  the  close  placement  of  the  fan  to  the  surface. 

Thermal  impedance  spectra  were  measured  on  fully- 
charged  Li-ion  batteries  manufactured  by  Panasonic,  Hita¬ 
chi,  Moly  Energy  and  Samsung  in  frequency  range  0. 1  Hz- 
1  mHz.  The  results  are  given  in  Fig.  6.  It  can  be  seen  that  the 
spectra  are  distinctly  different,  in  spite  of  the  similar  chem¬ 
istry  of  the  batteries.  This  observation  highlights  the  impor¬ 
tance  of  considering  thermal  properties,  characteristic  to 
each  particular  manufacturing  technology,  in  designing  heat 


management  of  large  battery  packs  with  high-power  cap¬ 
abilities. 

The  obtained  part  of  the  low-frequency  semicircle  was 
sufficient  to  determine  the  thermal  capacity  of  the  tested 
batteries.  In  order  to  determine  the  thermal  parameters,  the 
spectra  have  been  fitted  to  the  thermal  impedance  function  in 
Eq.  (5).  The  resulting  values  of  the  thermal  equivalent  circuit 
elements  are  given  in  Table  2  and  specific  values  are  in  listed 
Table  3. 

An  interesting  difference  is  observed  for  Moly  and  Sam¬ 
sung  batteries  compared  with  Panasonic  and  Hitachi  bat¬ 
teries.  The  latter  two  batteries  have  noticeable  thermal 
resistance  Rc  between  the  casing  and  the  electrode  stack, 
while  in  the  case  of  first  two  batteries  this  resistance  is 
negligible.  Absence  of  close  mechanical  contact  between  the 
electrode  stack  and  the  casing  could  be  reasons  for  this 
additional  thermal  resistance.  Such  loose  contact  should 
obviously  be  avoided  to  improve  the  protection  from  battery 
overheating.  The  lowest  specific  thermal  capacitance  is 
observed  for  the  Hitachi  battery,  which  also  has  the  lowest 
electric  capacity  as  shown  in  Table  1.  The  increase  of 


Fig.  6.  Comparison  of  thermal  impedance  spectra  of  Panasonic,  Moly  Energy,  Hitachi  and  Samsung  Li-ion  batteries  measured  between  0.1  Hz  and  1  mHz. 
Fit  lines  correspond  to  thermal  model  described  by  Eq.  (5)  with  parameter  values  given  in  Table  2. 
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Fig.  7.  Three-dimensional  plot  representing  dependence  of  thermal 
impedance  spectra  of  Panasonic  battery  on  state-of-charge.  Impedance 
spectra  are  obtained  in  frequency  range  from  0.1  Hz  to  1  mHz.  State-of- 
charge  ranges  from  discharged  state  (0%)  to  fully-charged  state  (100%). 


electric  capacity  in  battery  design  is  probably  related  to 
adding  more  components  with  high  thermal  capacity.  The 
electrode  stack  of  the  Hitachi  battery  also  has  the  highest 
thermal  conductivity  among  all  measured  batteries.  These 
two  properties  could  be  due  to  a  relatively  larger  contribu¬ 
tion  of  electrolyte  to  the  total  weight  of  battery. 

In  order  to  determine  how  sensitive  thermal  properties 
of  a  Li-ion  battery  are  with  respect  to  the  state-of-charge, 
thermal  impedance  spectra  of  Panasonic  Li-ion  battery 
were  obtained  at  different  states-of-charge.  The  battery 
was  fully  charged  and  then  subjected  to  a  repetitive  sequence 
of  discharge  at  the  5  h  rate  for  1  h,  relaxation  for  30  min 
and  impedance  spectrum  measurement  until  it  was  fully 
discharged.  The  impedance  spectra  obtained  at  different 
states  of  charge  is  shown  in  Fig.  7.  It  can  be  seen  that  the 
thermal  properties  of  Li-ion  battery  are  practically  indepen¬ 
dent  of  the  state-of-charge.  The  slight  increase  of  the  size  of 
the  semicircle  in  the  high-frequency  region  with  increasing 
state-of-charge  indicates  a  reduction  in  the  average  thermal 
conductivity  of  the  electrode  stack.  The  thermal  capacitance 
of  the  battery  remains  unchanged  during  discharge,  as  can  be 
seen  from  the  unchanged  semicircle  in  the  low-frequency 
region  of  the  thermal  impedance  plot. 


6.  Conclusions 

In  this  work,  we  have  demonstrated  a  possibility  of 
simultaneous  measurement  of  thermal  capacity  and  thermal 
conductivity  of  cylindrical  batteries  by  analysis  of  a  thermal 
pulse  response.  Derivation  of  the  thermal  impedance  spec¬ 
trum  from  the  temperature  response  transient  by  carrier 


function  Laplace  transform  has  been  described.  Thermal 
impedance  measurements  have  been  performed  on  Li-ion 
batteries  produced  by  Panasonic,  Moly  Energy,  Hitachi  and 
Samsung.  Considerable  differences  in  the  thermal  properties 
are  observed  despite  the  similar  basic  chemistry  and  design 
of  the  batteries. 

The  frequency-domain  thermal  impedance  function  for  a 
cylindrical  battery  have  been  derived  and  applied  to  the 
analysis  of  thermal  impedance  spectra.  Specific  thermal 
conductivities  and  capacities  have  been  calculated  for  all 
measured  batteries  by  fitting  experimental  thermal  impe¬ 
dance  spectra  to  the  theoretical  function.  These  values  can  be 
used  further  in  thermal  modeling  of  batteries,  which  is 
applied  for  design  of  high-power  systems  such  as  electric 
vehicles.  A  large  thermal  resistance  between  the  casing  and 
the  electrode  stack  detected  for  the  Panasonic  and  Hitachi 
batteries.  This  indicates  that  closer  fitting  of  the  electrode 
stack  can  improve  the  high-power  capabilities  of  these 
batteries. 

The  dependence  of  the  thermal  properties  on  the  state-of- 
charge  has  been  studied  for  a  Panasonic  Li-ion  battery. 
It  is  found  that  thermal  impedance  spectrum  is  virtually 
unchanged  during  battery  discharge.  A  slight  increase  in 
thermal  conductivity  with  state-of-charge  is  observed. 
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